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Abstract 

We consider the root lattice An and derive explicit formulae for the moments of its Voronoi cell. We then show 
that these formulae enable accurate prediction of the error probability of lattice codes constructed from j4„. 
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The root lattices An, Dn, Eq, Ej and Eg have attracted particular attention as structured codes for the AWGN 



^ channel yj. The highly symmetric structure of these lattices provides the grounds for extremely efficient decoding 
algorithms In this paper we consider codes constructed from the root lattice An and derive formulae for 

00 

accurately predicting the performance of these codes. This is achieved by deriving formulae for the moments of 
the Voronoi cell of An- Conway and Sloane suggested this approach to compute the quantizing constants (second 

j- ^ ' order moments) of the root lattices, |5]. In this paper we extend their technique to compute the moments of any 

HH order for An- 

^ In two dimensions A^ is the hexagonal lattice and in three dimensions is the face-centered cubic lattice. 
These are the densest known sphere packings in dimensions two and three and our results automatically include 
low dimensional codes constructed using these packings. In general, the lattice An does not produce asymptotically 
good codes in large dimensions, but does offer a coding gain in small dimensions. For these cases, we provide an 
error probability expression that can be computed to any degree of accuracy at any finite signal-to-noise ratio. 

This paper is organised as follows. In Section |ll] we give a brief overview of lattices and codes constructed from 
them, i.e., lattice codes. We describe lattice decoding and show how the probability of coding error can be expressed 
in terms of the moments of the Voronoi cell of the lattice used. Section |lll] states the main result (without proof), 
^ which gives recursive formula to compute the moments of the Voronoi cell of An- Section HVl describes the lattice 
^ An and some of its properties. An important property for our purposes is that the Voronoi cell of An is precisely 
^ the projection of a n + 1-dimensional hypercube orthogonal to one of its vertices [3^^. In Section IVl we use this 
K> property to show how integrals over the Voronoi cell of An can be expressed as integrals over the n-dimensional 
^ hypercube. These integrals are solvable and we use them to obtain the moments of the Voronoi cell in Section |Vll 
^ \ In Section I VII I we plot the probability of error versus signal to noise ratio for codes constructed from the lattices 
A\ ~ Z, A2, A3 ~ 1)3 and A4. We also plot the results of Monte-Carlo simulations that support our analytical 
results. 



II. Lattices, lattice codes, and lattice decoding 
A lattice. A, is a discrete subset of points in such that 

A = {x = Bu I u G Z"} 

where B € R'"^"^ is an m x n matrix of rank n, called the generator matrix or basis matrix or simply generator or 
basis. In particular, the set of n-tuples of integers Z" is a lattice (with the identity matrix as the generator) and we 
call this the integer lattice. A lattice A associated with a rank-n generator matrix B is said to be n-dimensional. 
If the generator is square, i.e. m = n, then the lattice points span and we say that the lattice is full rank. If 
B has more rows than columns, i.e. m > n, then the lattice points lie in a n-dimensional subspace of R"^. For 
any lattice A with an m x n generator matrix, we define 5a to be the hyperplane spanned by the columns of the 
generator matrix. It is easy to see that Sa is then the n-dimensional subspace containing the lattice A. 

The (open) Voronoi cell, denoted Vor(A), of a lattice A is the subset of Sa containing all points nearer (in 
Euclidean distance) to the lattice point at the origin than any other lattice point. It can be shown that the Voronoi 
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cell is an n-dimensional convex polytope that is symmetric about the origin. It is convenient to modify this definition 
of the Voronoi cell slightly so that the union of translated Voronoi cells UxeA Vor(A) + x is equal to 5a. That 
is, the Voronoi cell tessellates space when translated by points in A. To ensure this we require that if a face of 
Vor(A) is open, then its opposing face is closed. Specifically, if x € Vor(A) is on the boundary of Vor(A) then 
—X ^ Vor(A). We wont specifically define which opposing face is open and which is closed as the results that 
follow hold for any choice of open and closed opposing faces. 

The Voronoi cell encodes many interesting lattice properties such as the packing radius, covering radius, kissing 
number, minimal vectors, center density, thickness, and the normalized second moment (or quantizing constant) 111 
0]. The error probability of a lattice code can also be evaluated from the Voronoi cell as we will see. There exist 
algorithms to completely enumerate the Voronoi cell of an arbitrary lattice iF/l-floll. In general these algorithms are 
only computationally feasible when the dimension is small (approximately n < 9). Even with a complete description 
of the Voronoi cell it is not necessarily easy to compute the probability of coding error. 

The Voronoi cell is linked with the problem of lattice decoding. Given somepoint y S M" a lattice decoder (or 
nearest lattice point algorithm) returns the lattice point in A that is nearest to y llllll . Equivalently it returns the lattice 
point X such that the translated Voronoi cell Vor(A) +x contains y. Computationally lattice decoding is known to 
be NP-hard under certain conditions when the lattice itself, or rather a basis thereof, is considered as an additional 



input parameter lll2l Il3ll . Nevertheless, algorithms exist that can compute the nearest lattice point in reasonable 



time if the dimension is small (approximately n < 60). One such algorithm is the sphere decoder jllL Il4l Il5ll . 



A good overview of these techniques is given by Agrell et. al. Illll . Fast nearest point algorithms are known for 
specific lattices mi Il6ll . For example, the root lattices and An&nd their dual lattices and ^* can be 



decoded in linear-time, i.e. in a number of operations of order 0(n) yl, |3[]. 

Lattices can be used to construct signal space codes. A signal space code X of dimension n is a finite set of vectors 
(or points) in the Euclidean space M". Each vector in X is called a codeword and represents a particular signal. 
The number of codewords is denoted by \X\. If each codeword is transmitted with equal probability then the rate 
of the code is i? = ^ log2 \X\ bits per codeword. The average power of the code is given by P = |^ SxeA" ll^lP- 
In the AWGN channel the received signal takes the form 

y = X + w 

where y G R", x G and w is a vector of independent and identically distributed Gaussian random variables with 
variance cr^. If the receiver employs maximum likelihood decoding then the estimator of x given y at the receiver 
is 

X = argmin ||y — x||^. (1) 

That is, the receiver computes the codeword in X nearest in Euclidean distance to the received signal y. Assuming 
that each codeword is transmitted with equal probability then the probability of correct decoding is 

and the probability of error is 

A lattice code is a signal space code with codewords taken from a finite subset of points of some lattice A in 
M". There are infinitely many ways to choose a finite subset from a lattice, but common approaches make use of 
a bounded subset of R", called a shaping region S C R". The codewords are given by those lattice points inside 
the shaping region, that is, 

x = Sr^K. 

Common choices of shaping region are n-dimensional spheres, spherical shells, hypercubes, or the Voronoi cell 
of a sublattice of A rl7l-ll9ll. A consequence of the finiteness of a lattice code is that the maximum likelihood 
decoder ([T]) is not equivalent to lattice decoding. Computing a nearest lattice point will not in general return a 
lattice point from the code (the decoded lattice point might Lie outside the shaping region). 
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Let X he a lattice code constructed from a lattice A. Assuming that each codeword x G has an equal probability 
of being transmitted and that the receiver employs lattice decoding, then the average probability of correct decoding 
in the AWGN channel is 

Pc < PcW = lim y Pr(x + w e Vor(A) + x) 

lA"!— s>oo \X\ ^ — ' 

= Pr(x € Vor(A)) 

= [ e-ll'^IIV^'^^dx. (2) 

(V2vro-)" Jvor(A) 

The probability of error is Pe < 1 — i-'c'(A), where the upper bound is asymptotically tight for large constellations 
(the proportion of the codewords near the boundary of the shaping region becomes small). By expanding = 
1 + x + ^ + . . . according to its Maclaurin series we obtain 



^— T ^-'^^ f llxll-dx 



(3) 



So, to obtain arbitrarily accurate approximations to the probability of error it is enough to know the values of 
/vor(A) l|x|P™(ix for m = 1, 2 . . . for some sufficiently large m. The number of terms required will increase as the 
noise variance gets smaller. This implies that the bound with a fixed number of terms is not asymptotically tight 
but is very accurate up to a finite signal-to-noise ratio. 

In this paper we focus on n-dimensional lattice codes constructed from the family of lattices called An and we 
will derive expressions for the mth moments 



Mn{m) = [ ||xf '"dx. 



/Vor{A„) 

These can be summed in ^ to give arbitrarily accurate approximations for the probabiUty of error. 

III. The main result 
We now state our main result. The mth moment M„(m) of the lattice An satisfies 



Mn{m) n^n + 1 y-? G(n — 1, a, 2k — 2a — b) 

ml n + 2m ^^^^^ {—1)^~"-H{n,m,k,a,b) 



k=0 a=0 6=0 

where the function 

(n + l)"'-''a\(m- ky.b\(k - a-bV. 

H[n,m,k,a,b) = r ; . 

^ ' ' ' ' y 2"fi"^~''^ 

and the function G{n, c, d) satisfies the recursion 



f c\f d\ G{n -l,c-c',d- d') 

G{n, c,d) = Y^J2 J [d> J ^TTdTTi 

c'=Od'=0 ^ / ^ / 



with the initial conditions 

G(l, c, d) = and G(n, 0, 0) = 1. 

^ ' ^ 2c + d+l V > > ; 
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For fixed m it is possible to solve this recursion in n and obtain formula for the Mn{m) in terms of n. The first 
five such formula are: 



M„(0) = Vn + 1 the volume of Vor(A„), 
n(n + 3) 

, the second i 

50n + 55n2 + 34n3 + 5n^ 



M„(l) = the second moment of Vor^„ [Jj, p. 462], 

12vn + 1 



M„(2) 
M„(3) 
M„(4) 



720(1 + n)3/2 
1960n + 2142^2 + 2681n3 + 1423n^ + 399n^ + 35n6 
60480(1 + n)5/2 ' 
93744n + 34356^2 + 112172^3 + 89343n'^ + 5322471^ + 1724671*^ + 2940n^ + 17571^ 



3628800(1 + n)7/2 

We have explicitly tabulated these formula for m = to 40. For larger m direct evaluation for specific n from the 
recursive formula is preferable. We will derive these results in Section |Vll but first need some properties of the 
lattice An- 

IV. The lattice An 

Let H be the hyperplane orthogonal to the all ones vector of length n + 1, denoted by 1, that is 

1=[1 1 ••• 1]', 

where superscript ' indicates the transpose. Any vector in H has the property that the sum (and therefore the mean) 
of its elements is zero and for this reason H is often referred to as the zero-sum plane or the zero-mean plane. 
The lattice An is the intersection of the integer lattice with the zero-sum plane, that is 

yl„ = Z"+ini/ = {xGZ"+i |x'l = 0}. (6) 

Equivalently, An consists of all of those points in Z"+^ with coordinate sum equal to zero. The lattice has n(n + 1) 
minimal vectors, each of squared Euclidean length 2, so the packing radius is The n-volume of the Voronoi 
cell Vor(A„) is y/TT+T |1, p. 108]. 

The Voronoi cell of An is closely related to the n + 1 dimensional hypercube Vor(Z"+^) as the next theorem 
will show. This result has appeared previously Jsl @], but we repeat it here so that this paper is self contained. We 
denote by 

Q = - - 



I'l n + 1 

the projection matrix orthogonal to 1 (i.e. into the zero-sum plane) where I is the n + 1 by n + 1 identity matrix. 
Given a set S of vectors from R"+^ we write QS to denote the set with elements Qs for all s € S*, i.e. the set 
containing the projection of the vectors from S. 

Lemma 1. The projection o/ Vor(Z"+^) into the zero-sum plane is a subset o/ Vor(A„). That is, 

QVor(Z"+i) C Vor(A„). 

Proof: Let y G Vor(Z"+^). Decompose y into orthogonal components so that y = Qy + 1\ for some f G R. 
Then Qy G Q Vor(Z"+^). Assume that Qy ^ Vor(A„). Then there exists some x G An such that 

l|x-Qy|P < ||0-Qy|p ^ ||x-y + tl|p < ||y - tlf 

=> ||x - y|p + 2tx'l < ||y|p. 

By definition © x'l = so ||x - y|p < ||y|p. This violates that y G Vor(Z"+^) and hence Qy G Vor(A„)Q ■ 

'This proof can be generalised to sliow tliat for any lattice L and hyperplane P such that PnL is also a lattice it is true that pVor(L) C 
Vor(L n P) where p indicates the orthogonal projection into P Lemma 2.1]. 
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Theorem 1. The projection o/'Vor(Z"+^) into the zero-sum plane is equal to Vor(A„). That is, 

Vor(A„) = QVor(Z"+i). 

Proof: Let denote a vector with ith element equal to one and the remaining elements zero. The n-volume of 
Vor(A„) is \Jn + 1. From Berger et. al. [20, Theorem 1.1] we find that the n-volume of the projected hypercube 

QVor(Z"+i) is equal to 

n+l , n.+ l 

1=1 1=1 ^ 

also. It follows from Lemma [T] that Q Vor(Z"+^) C Vor(A„), so, because the volumes are the same, and because 
Vor(A„) and QVor(Z"+i) are polytopes, we have Vor(A„) = QVor(Z"+i). ■ 

V. Integrating a function over Vor(A„) 

We would like to be able to integrate functions over the Voronoi cell of An- Consider a function / : M"+^ i— )• M. 
The definition we have made for An above places it in the n-dimensional zero-sum plane, lying in M"+^. The Voronoi 
cell is a subset of the zero-sum plane that has zero n + 1-dimensional volume. So, the integral J^orlA ) /(^)^^ 
equal to zero. This is not what we intend. By an appropriate change of variables it would be possil5le to write the 
Voronoi cell Vor(j4„) in an n-dimensional coordinate system, and then integrate. However, we find the following 
approach simpler. Given a set S of vectors from the zero sum plane, let 5 x 1 denote the set of elements that can 
be written as x + y where x € S" and y = kl for some k G [—1/27 V^]- Now, the integral over the Voronoi cell can 
be written as 



/(Qx)dx= / /(Qx)dx. (7) 

Vor(yl„)xl JQVor{Z"+i)xl 

It is not immediately clear how an integral over Vor(j4„) x 1 should be performed. Consider the following simpler 
integral over the hypercube Vor(Z"+^), 

/ /(Qx)dx. (8) 

JVor(Z"+i) 

This integral is not equal to (|7]) because, although Qx is always an element of Vor(j4„), the integral is not uniform 
over Vor(A„). To see this, consider some x G Vor(Z"+^) and let Xmax be the maximum element of x and be 
the minimum element. Then x + A;l G Vor(Z""'"^) for those k G [— 1/2 — Xmin, ^1'^ — Xmm)- The length of this interval 
is 1 — Xmax + SO the (onc dimensional) volume of points in Vor(Z"+^) that, once projected orthogonally to 1, 
are equal to Qx is given by 



min ) 



II 1 11(1 - 2;max + a^min) = + 1(1 - X^ax + 

The integral ([7]) can be obtained by normalising ([8]l by this length, that is, 

/ /(Qx)dx= / ^ ^ -^^Q""^ -dx. 

iVor(A„)xl JVor(Z"+i) \/n + 1(1 - Xmax + a^min) 

The primary advantage of this integral is that the bounds are given by the hypercube Vor(Z"+^). 

Let us now restrict /(x) so that it depends only on the magnitude ||x||, for example /(x) = Hxp™ could be 
a power of the Euclidean norm of x. Now /(x) is invariant to permutation of x. Let x be such that xi is the 
maximum element and X2 is the minimum element. Our integral is now equal to 

n(n+i) c^i^ r r r /(q^) , , , 

-axn+i ■ ■ ■ dx2 dxi. 



Vn + 1 J -1/2 J -1/2 Jx2 Jx2 'i--Xl + X2 

The factor n(n + 1) arises because there are n(n + 1) ways to place two elements (i.e. xi and X2) into n + l 
positions. 

We can make further simplifications. Letting t = xi — X2 and y = xi + 1/2 and changing variables gives 

/■I ry /-y-i/s ry-^/'^ /(Qx) 

nVn + 1 / / •••/ —dxn+i.-.dx^dtdy, 

Jo Jo Jy-t-1/2 Jy-t-1/2 i — i 
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and letting 11)1-2 = — y + 1/2 + t for i = 3, . . . , n + 1 gives 

.1 i-y ft 



nVnTT [ r [ ■■■ [ ^^^^dwn-i ... dwidtdy. (9) 
JO Jo Jo Jo 1 — i 

Observe that x = w + (7/ — t — 1/2)1 where w is the column vector 

W = [t,0,Wl,W2, . . . ,Wn-l]'- 

Projecting orthogonal to 1 gives Qx = Qw. Interestingly w does not contain y so the term inside the integral does 
not depend on y. This is the integral we will use to compute the moments of An- 

Example 1. (The volume of the Voronoi cell) In order to demonstrate this approach we will derive the 0th moment 
(i.e. the volume) of the Voronoi cell using Setting /(Qw) = ||Qw||° = 1 we obtain, 

/■I ry f-t rt i 

Mn{0) = n\Jn + 1 / / / ■ ■ ■ / dwn-i ■ ■ ■ dwi dt dy 

Jo Jo Jo Jo ^ — t 



lO JO JO Jo 
cy fl fn-1 



Vn + 1 / / -dtdy. 

Jo Jo ~ * 

Vn + l / f3{y,n,0)dy = y/n + l. 
Jo 



as required. Here I3{x,a,b) = Jq t"- ^(1 — t)^ ^ dt is the incomplete beta function 112 ill and we have used the 
identity (3{y, n, 0) dy = i. 

VI. The moments of 
We now derive expressions for the Mn{m). Setting /(Qx) = (HQwlp)"" in Q we obtain, 

n^Jn + l Jo Jo Jo Jo 
Now IIQwp = ||w|p — ;^^(w'l)^ and recalling that w = [t,0,wi, . . . ,Wn-i]' we can write 

||Qw||2 = ||wf-^(w'l)2 

n + 1 



n+l ^ / n+1 

i=l \ 1=1 

n+l 



^ I n+l /n+l \2\ 



say, where 



and where. 



C = [ f and D = A- -^B - -^B^ 

n + ly n + l n + l 



n—l n—1 

Wi 

i=l 1=1 



A = '^wf and ^ = XI 

i=l 

Mn{m) fy /■* {C + D 



Now, 

nVn + 1 ^0 ^0 ^ Jo 1 - * 
and by expanding the binomial (C + D)"^ we get 

M„(m) fy 1 



dwn~i ■ ■ ■ dvji dt dy, 



n 
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ki+k2+k3=k 



Expanding as a trinomial gives 

^ k\ 

n2k-2a-b / i \ 

= V V ( ~ ] 2H^ 

^ ^ a\b](k -a -by. \n + 1 J 

where the second Une follows by setting ki = a, k2 = b and k^ = k — a — b. In Appendix |A| we show that the 
integral of A"' B'^^~'^"-~^ over wi, . . . Wn-i is 

/* • • • /* A^B^''^^^-^dwn-i ...dwi = e-^+^''-^G{n - 1, a, 2k -2a- b). (10) 
JO Jo 

where G(n, c, d) satisfies the recursion given by ([5]). So, let P satisfy 

P = fi-"-2'= /"*... /"* D^dwn-i ...dwi 
Jo Jo 

_j^j^2'>k\G{n-l,a,2k-2a-b)f -1 ^ '^-^ 



a=0 b=0 ^ ' ^ 



/ \ m—k 

Now C""-'' = (irh) t2(m-fe) ^j^j 



1 \ ^ / m\ / n 



n + 2m ^-^ \ k I \ n+l 

k=o ^ ^ ^ 



This expression is equivalent to that from (H]). 

VII. Results and simulations 

We now plot the probability of coding error versus signal to noise ratio (SNR) for the lattices Ai, A2, A'^, A^, 
and A^. For these plots the SNR is related to noise variance according to ^ 

y2/n 

where V is the volume of the Voronoi cell and n is the dimension of the lattice. Figure [T] shows the 'exact' 
probability of error (correct to 16 decimal places) computed using the moments Af„(m) and ([3]l (solid line). The 
number of moments needed to ensure a certain number of decimal places accuracy depends on n and also on 
the noise variance cr^. At most 321 moments where needed for Figured] We also display the probability of error 
computed approximately by Monte-Carlo simulation (dots). The simulations are iterated until 500 error events occur. 

The plot also display an approximation for the probability of error for the 8-dimensional E'g lattice. The 
approximation is made in the usual way by applying the union bound to the minimal vectors of the lattice iH 
p. 71]. The £^8 lattice has 240 minimal vectors of length \/2. The packing radius of £^8 is therefore p = \/2/2. 
Applying the union bound the probability of error satisfies 

P,< 240erfc(-|-)= 240erfc(± 

where erfc(x) = 1 — erf (x) is the complementary error function. For the lattice this approximation is an upper 
bound because the relevant vectors of (those vectors that define the Voronoi cell) are precisely the 240 minimal 
vectors. 
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Fig. 1. The probability of error versus SNR for Ai ~ Z, A2, A3 ~ D3, Ai, A5, Ag and -Eg. 



VIII. Conclusion 

Recursive formulae for the moments of the Voronoi cell of the lattice An were found. These enable accurate 
prediction of the performance of codes constructed from A„. In two dimensions A2 is the hexagonal lattice and in 
three dimensions ^3 is the face-centered cubic lattice. Our results include codes constructed using these packings 
as a special case. 

Appendix 

A. A multinomial type integral over a hypercube 
In ( fTOl ) we required to evaluate integrals of the form 

F{n-l,a,2k-2a-h)= [ ■■■ [ A'^B^'"'^'''^ dwn-i ■ ■ ■ dwi, 

Jo Jo 

or equivalently, integrals of the form 




where n, c and d are integers. We shall find a recursion describing this integral. Write 
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Expanding the two binomials gives 



EE/* -i'O)©^""^' ' 



c'=0(i'=0''" N / \ / yj=l y 

Integrating the Xn term gives 



- - . d-d' 



Note that 



F(n — 1, c — c', d — c?') = / ••• / y^a^ll 5^^* dxi ■ ■ ■ dxn-i- 

Jo Jo \ ,=1 / Vi=i / 



So F{n, c, d) satisfies the recursion 



Fin,.. = E E (^-) - - 

c'=Od'=0 \ / \ / 

with the initial conditions 

i2c+d+l 

The F(n, c, d) can be written as t"+^^+'^G(n, c, d) where G(n, c, d) is rational. To see this write 

'd\ t2c'+d'+l 



Fin, = E E (^,) - M - ^, « - ^) 

c'=Od'=0 \ / \ / 

^ \d' y 2c' + + 1 



c'=Qd'=Q 

Now G(n, c, d) is the rational number satisfying the recursion 



g(„.c,.)^ee(:)(:,) 



c\ ( d\ G{n — 1, c — c' , d — d') 



c'=0 d'=0 

with the initial conditions 



2c' + d' + l 



6. Solving this recursion for fixed d and c 

For fixed d and c this recursion can be solved explicitly. Write 

'd\G{n-l,c-c',d-d') 



G{n,c,d) = G{n-l,c,d)+ ^ fulfil 



2c' + d' + 1 



(c',<i')^(0,0) 

where the sum Yl{c' d')^(o o) ^^^^ all < c' < c and < d' < d except when both d and c are zero. Denote by 
Q{z,c,d) the z-transform of G{n,c,d). Taking the z-transform of both sides in the equation above gives 

sr-^ f c\ f d\Q{z,c — c' ,d — d') 



(c',d')^(0,0) 



2c' + d' + 1 
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So the z-transform Q{z, c, d) satisfies this recursive equation. The initial condition is G{z, 0, 0) = jr^- By inverting 
this z-transform and using the resultant expressions in ^ we obtain formulae in n for the moment Mn{m). This 
procedure was used to generate the formula described in Section JII] Mathematica 8.0 was used to perform these 
calculations. 
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